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Abstract 

For an nt transmit, rir receive antenna system (nt x rir system), a full-rate space time block code 
(STBC) transmits n,„i„ = ■min{nt, rir) complex symbols per channel use. The well known Golden code 
is an example of a full-rate, full-diversity STBC for 2 transmit antennas. Its ML-decoding complexity is 
of the order of Af for square A/-QAM. The Silver code for 2 transmit antennas has all the desirable 
properties of the Golden code except its coding gain, but offers lower ML-decoding complexity of the 
order of A/^. Importantly, the slight loss in coding gain is negligible compared to the advantage it offers in 
terms of lowering the ML-decoding complexity. For higher number of transmit antennas, the best known 
codes are the Perfect codes, which are full-rate, full-diversity, information lossless codes (for rir > rit) 
but have a high ML-decoding complexity of the order of A/"'"™'" (for rir < rit, the punctured Perfect 
codes are considered). In this papej], a scheme to obtain full-rate STBCs for 2° transmit antennas and 
any rir with reduced ML-decoding complexity of the order of is presented. The codes 

constructed are also information lossless for rir > rit, like the Perfect codes and allow higher mutual 
information than the comparable punctured Perfect codes for rir < nt - These codes are referred to as the 
generalized Silver codes, since they enjoy the same desirable properties as the comparable Perfect codes 
(except possibly the coding gain) with lower ML-decoding complexity, analogous to the Silver-Golden 
codes for 2 transmit antennas. Simulation results of the symbol error rates for 4 and 8 transmit antennas 
show that the generalized Silver codes match the punctured Perfect codes in error performance while 
offering lower ML-decoding complexity. 

Index Terms 

Anticommuting matrices, ergodic capacity, full-rate space-time block codes, low ML-decoding com- 
plexity, information losslessness. 

'Part of the content of this manuscript has been presented at IEEE ISIT 2010 and another part at IEEE Globecom, 2010. 
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I. Introduction and Background 

Complex orthogonal designs (CODs) [1], [2], although provide linear Maximum Likelihood (ML)- 
decoding, do not offer a high rate of transmission. A full-rate code for an nt x MIMO system transmits 
min{nt,nr) independent complex symbols per channel use. Among the CODs, only the Alamouti code 
for 2 transmit antennas is full-rate for a 2 x 1 MIMO system. A full-rate STBC can efficiently utilize 
all the degrees of freedom the channel provides. In general, an increase in the rate tends to result in an 
increase in the ML-decoding complexity. The Golden code [3] for 2 transmit antennas is an example of 
a full-rate STBC for any number of receive antennas. Until recently, the ML-decoding complexity of the 
Golden code was reported to be of the order of M^, where M is the size of the signal constellation. 
However, it was shown in [4], [5] that the Golden code has a decoding complexity of the order of A'P-^ 
for square M-QAM. Current research focuses on obtaining high rate codes with reduced ML-decoding 
complexity (refer to Sec. JI] for a formal definition). For 2 transmit antennas, the Silver code, named so 
in [6], was first mentioned in [7] and independently presented in [8] along with a study of its low ML- 
decoding complexity property. It is a full-rate code with full-diversity and an ML-decoding complexity 
of the order of for square Af-QAM. Its algebraic properties have been studied in [6] and [9] and a 
fixed point fast decoding scheme has been given in [10]. For 4 transmit antennas, Biglieri et. al. proposed 
a rate-2 STBC which has an ML-decoding complexity of the order of M^-^ for square Af-QAM without 
full-diversity [11]. It was, however, shown that there was no significant reduction in error performance 
at low to medium SNR when compared with the previously best known code - the DjABBA code [12]. 
This code was obtained by multiplexing Quasi-orthogonal designs (QOD) for 4 transmit antennas [13]. 
In [4], a new full-rate STBC for 4 x 2 system with an ML-decoding complexity of M^'^ was proposed 
and was conjectured to have the non-vanishing determinant (NVD) property. This code was obtained by 
multiplexing the coordinate interleaved orthogonal designs (CIODs) for 4 transmit antennas [14]. These 
results show that codes obtained by multiplexing low complexity STBCs can result in high rate STBCs 
with reduced ML-decoding complexity and by choosing a suitable constellation, there won't be any 
significant degradation in the error performance when compared with the best existing STBCs. Such an 
approach has also been adopted in [15] to obtain high rate codej^ from multiplexed orthogonal designs. 
More recently, full-rate STBCs with an ML-decoding complexity of the order of M^-^ and a provable 
NVD property for the 4x2 system have been proposed in [19] and [20]. 

'Fast decodable STBCs have been constructed in [16]-[18], but these codes are not full-rate in general, and make use of near 
ML-decoding algorithms. 
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In general, it is not known how one can design full-rate STBCs for an arbitrary number of transmit 
and receive antennas with reduced ML-decoding complexity. It is well known that the maximum mutual 
information achievable with an STBC is at best equal to the ergodic capacity of the MIMO channel, 
in which case the STBC is said to be information lossless (see Section |ll] for a formal definition). It 
is known how to design information lossless codes [21] for the case where > nj. However, when 
rir < nt the only known code in literature which is information lossless is the Alamouti code, which 
is information lossless for the 2x1 system alone. It has been shown in [7], [22] and [12] that when 
rir < nt, self-interference of the STBC (a formal definition of self interference is given in Section has 
to be minimized for maximizing the mutual information achieved with the STBC. Not much researcq; 
has been done on designing codes that allow a high mutual information when < nt. In this paper, for 
nt = 2", we systematically design full-rate STBCs which have the least possible self-interference and 
the lowest ML-decoding complexity among known full-rate STBCs for < nt and consequently allow 
higher mutual information than the best existing codes (the Perfect codes with puncturing [25], [26]), 
while for > nt, the proposed STBCs are information lossless like the comparable Perfect codes. We 
call these codes the generalized Silver codes since, analogous to the silver code and the Golden code 
for 2 transmit antennas, the proposed codes have every desirable property that the Perfect codes have, 
except the coding gain, but importantly, have lower ML-decoding complexity. The contributions of the 
paper are: 

1) We give a scheme to obtain rate-1, 4-group decodable codes (refer Section HI] for a formal definition 
of multi-group decodable codes) for nt = 2" through algebraic methods. The speciality of the 
obtained design is that it is amenable for extension to higher number of receive antennas, resulting 
in full-rate codes with reduced ML-decoding complexity for any number of receive antennas, unlike 
the previous constructions [27]-[29] of rate-1, 4-group decodable codes. 

2) Using the rate-1, 4-group decodable codes thus constructed, we propose a scheme to obtain the 
generalized Silver codes, which are full-rate codes with reduced ML-decoding complexity for 2" 
transmit antennas and any number of receive antennas. These codes also have the least self- 
interference among known comparable STBCs and allow higher mutual information with lower 
ML-decoding complexity than the comparable punctured Perfect codes for the case 7ir < nt, 

^The full-rate STBCs in [23], designed for Ur < nt, are not linear dispersion codes. They are based on maximal orders and 
use spherical shaping due to which the encoding and decoding complexity is extremely high. The STBCs in [24], also designed 
for Ur < nt, use the concept of restricting the number of active transmit antennas to be no larger than the number of receive 
antennas, and so, the mutual information analysis for these codes is very difficult. These STBCs are diversity-multiplexing gain 
tradeoff (DMT) optimal but are associated with a very high ML-decoding complexity. 
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while being information lossless for > nt. In terms of error performance, by choosing the 
signal constellation carefully, the proposed codes have more or less the same performance as the 
corresponding punctured Perfect codes. This is shown through simulation results for 4 and 8 transmit 
antenna systems. 

The paper is organized as follows. In Section |II1 we present the system model and the relevant 
definitions. The criteria for maximizing the mutual information with space time modulation are presented 
in Section JII] and our method to construct rate- 1 , 4-group decodable codes is proposed in Section |IVl 
The scheme to extend these codes to obtain the generalized Silver codes for higher number of receive 
antennas is presented in Section |Vl Simulation results are discussed in Section |Vl] and the concluding 
remarks are made in Section IVIII 

Notations: Throughout, bold, lowercase letters are used to denote vectors and bold, uppercase letters 
are used to denote matrices. Let X be a complex matrix. Then, and X^ denote the Hermitian and 
the transpose of X, respectively and unless used to denote indices or subscripts, j represents The 
entry of X is denoted by X(i, j) while tr(X) and (iet(X) denote the ti^ace and determinant of X, 
respectively. The set of all real and complex numbers are denoted by M and C, respectively. The real 
and the imaginary part of a complex number x are denoted by x/ and xq, respectively. ||X|| denotes 
the Frobenius norm of X, ||x|| denotes the vector norm of a vector x, and It and Ot denote the T x T 
identity matrix and the null matrix, respectively. The Kronecker product is denoted by ^ and vec(X.) 
denotes the concatenation of the columns of X one below the other. For a complex random variable 
X, K[X] denotes the mean of X and Kx {f{X)) denotes the mean of f{X), a function of the random 
variable X. The inner product of two vectors x and y is denoted by (x, y). For a set S, aS = {as\s G S}. 
Let V and Q be two sets such that P D Q. Then V\Q denotes the set of elements of V excluding the 
elements of Q. For a complex variable x, the (.) operator acting on x is defined as 



^ A 
X = 



XI -XQ 
XQ XI 



The (.) can similarly be applied to any matrix X G c^xm by replacing each entry Xij with Xij, i = 
1,2,- •• ,n,j = 1,2,- •• ,m, resulting in a matrix denoted by X G ]^2nx2m Qjven a complex vector 
X = [xi,X2,--- ,Xn]'^, X is defined as x = [xii,xiq, ■ ■ ■ ,Xni,XnQ]'^ ■ It follows that for A G ([^mxn^ 
B G C"^P and C = AB, the equalities C = AB and i^^^C) = (Ip ® A)u1^) hold. 
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II. System Model 

We consider the Rayleigh block fading MIMO channel with full channel state information (CSI) at 
the receiver but not at the transmitter. For nt x MIMO transmission, we have 



where S G c^txT jj^g codeword matrix whose average energy is given by E(||S|p) = ritT, N G C"--^^ 
is a complex white Gaussian noise matrix with i.i.d. entries ~ Mc (0, 1) (complex normal distribution 
with zero mean and unit variance), H G xnt jj^g channel matrix with the entries assumed to be 
i.i.d. circularly symmetric Gaussian random variables ~ J\fc (0,1), Y G (^rirxT jj^g received matrix 
and SNR is the signal-to-noise ratio at each receive antenna. 

Definition 1: {Code rate) Code rate is the average number of independent information symbols trans- 
mitted per channel use. If there are k independent complex information symbols (or 2k real information 
symbols) in the codeword which are transmitted over T channel uses, then, the code rate is k/T complex 
symbols per channel use {2k /T real symbols per channel use). 

Definition 2: {Full-rate STBCs) For an rit x n^, MIMO system, if the code rate is min {nt,nr) complex 
symbols per channel use, then the STBC is said to be fiill-rate. 

Assuming ML-decoding, the metric that is to be minimized over all possible values of codewords S 
is given by 



Definition 3: {ML-Decoding complexity) The ML decoding complexity is measured in terms of the 
maximum number of symbols that need to be jointly decoded in minimizing the ML decoding metric. 
For example, if the codeword transmits k independent symbols of which a maximum of p symbols need 
to be jointly decoded, the ML-decoding complexity is of the order of AP, where M is the size of the 
signal constellation. If the code has an ML-decoding complexity of order less than M'^, the code is said 
to have reduced ML-decoding complexity. 

Definition 4: {Generator matrix) For any STBC that encodes 2k real symbols (or k complex informa- 
tion symbols), the generator matrix G G K2Tntx2fc defined by [11] 




(1) 




vec (S) = Gs, 



where S is the codeword matrix, s = [si, S2, • " ' 



S2k]^ is the real information symbol vector. 
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A codeword matrix of an STBC can be expressed in terms of weight matrices (linear dispersion 
matrices) [30] as 

2k 

S = ^ ^ SjAj. 

i=l 

Here, Ai,i = 1,2,- •• ,2k, are the complex weight matrices of the STBC and should form a linearly 
independent set over M. It follows that 

vec{Ai) vec{A2) ■■■ vec{A2k] 

Due to the constraint that E(||Sp) = n^T, we have, X]^=i IE(si)^tr (AjAf ) = mT. Choosing 
E{sif = 1/2 for alH = 1, • • • , 2k, we have 

2k 

Y,tr{AiAf)=2ntT. (2) 

1=1 

Definition 5: {Multi-group decodable STBCs) An STBC is said to be (7-group decodable [29] if its 
weight matrices can be separated into g groups Qi, Q2, • • • , Gg such that 

A,Af +AjAf = 0„,, AiGGi, AjGGp, l,p £ {1,2,- ■ ■ ,g}, I ^ p. 
Definition 6: {Self-interference) For an STBC given by S = Yl^=i ^i^i^ the self-interference matrix 
[12] is defined as 

2k-l 2k 
i=l j>i 

Definition 7: {Punctured Codes) Punctured STBCs are the codes with some of the symbols being 
zeros, in order to meet the full-rate criterion. 

For example, a codeword of the Perfect code for 4 transmit antennas [25] transmits sixteen complex 
symbols in four channel uses and has a rate of 4 complex symbols per channel use. If this code were 
to be used for a two receive antenna system which can only support a rate of two independent complex 
symbols per channel use, then, eight symbols of the Perfect code can be made zeros, so that the codeword 
transmits eight complex symbols in four channel uses. These eight symbols correspond to the two layers 
[25] of the Perfect code. 

Equation ([T]) can be rewritten as 



' — ■ / SNR ' — 

vec(Y) = \ / He„s + vecCN) , (3) 

V nt 

where Hgg £ ]^2nrTx2n^^„T ^ called the equivalent channel matrix, is given by Heq = (It "X" H) G, with 
G e ]^2ntTx2n,„,„T bgjj^g (^e generator matrix as in Definition IH 
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Definition 8: {Ergodic capacity) The ergodic capacity of an nt x n^. MIMO channel is [31] 



C„,xn. = Eh {log det {in^ + ^^HH^^^ 



With the use of an STBC, the maximum mutual information achievable is [32] 

'^STBC = {log det {l2nrT + '^^"^ HegHf^^ ^ . 

It is known that C^^xn^ > '^stbc- If Cruxn^ = ^stbc^ the STBC is said to be information lossless. If 
the generator matrix G is orthogonal (from Definition IH this case arises only if > nt and the STBC 
is full-rate, i.e, k = ntT), the STBC is information lossless. 

III. Relationship between weight matrices and the maximum mutual information 

Capacity can be achieved with the use of continuous inputs with Gaussian distribution. If one were 
able to use continuous Gaussian distributed inputs in practice, using the V-blast scheme would suffice, 
since diversity is irrelevant. But in practice, one has to use finite discrete inputs, and diversity becomes an 
important aspect, necessitating the use of full-diversity STBCs. Even though we considered the limited 
block length scenario for space-time coding as a standalone scheme, in practice, one would also have an 
outer code and coding would be done over large block lengths to go close to capacity. In such a scenario, 
the maximum mutual information that an STBC allows becomes an important parameter for the design of 
STBCs. It is preferable to use STBCs which allow mutual information as close to the channel capacity as 
possible. It has been shown that if the generator matrix is orthogonal, the maximum mutual information 
achievable with the STBC is the same as the ergodic capacity of the MIMO channel [21], [32]. For the 
generator matrix to be orthogonal, a prerequisite is that the number of receive antennas should be at least 
equal to the number of transmit antennas. When < nt, only the Alamouti code has been known to 
be information lossless for the 2x1 MIMO channel. In [22], by using the well known matrix identities 
det M = e*'^('°^''^) and log{l + X) = Y.n=i ^"^^"'' x", an expansion of the ergodic MIMO capacity in 
SNR was obtained as 

oo 
n=l 

with C„ = ^ (^) \.^^^)^Y) ■ "^^^ ^^^^ "-^^ coefficients can easily be checked to be Ci = 

and C2 = —nr{nr + nt)/nt. On a similar note, Istbc can also be expanded in SNR as Istbc = 
En=i InSNR', where 
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(HegH 



— Eh Ur 



(H^gHgg) 



(4) 



2Tn \ rit 

Let H = H^gHgg. It is straightforward to check that H(z, j) = itr (S^H^H), where S^^- = AjA^^ + 
AjA^ . Hence, 

2Tn 



h = ^EH(tr[Hr,H,,]) = ^5]^EH(tr(S..H^H)) 

i=l 

^5]tr (A,AfEH (H^H))=n„ 



where E (H^H) = rirlnt and ^ is used in obtaining Ii. So, using all the available power helps one to 
achieve the first order capacity. The second coefficient I2 has been calculated in [7] to be 



J2 E (iKS|) + n.(trS.,)2). (5) 



IGTnf 

i=l j=i 

In [7], it was argued that typical discrete input schemes fail to achieve capacity at the third order in 
the expansion of the mutual information and hence, I2 should be maximized. From ([S]), it is clear that 
to maximize I2, the following criteria should be satisfied. 

1) Hurwitz-Radon Orthogonality: as many of Sjj should be equal to O^^ as possible, for 1 < i < j < 
2Tnr. 

2) Tracelessness: Sij should be traceless, for all 1 < i < j < 2Tnj.. 

In fact, the first criterion, which is equivalent to minimizing the self-interference, is already clear from 
(111), where it can be observed that a larger number of zero entries of HjgHgq contributes to a lower value 
of the trace of (H^gHgg)^. Hence, to design a good STBC with a high mutual information when < nt, 
one should have as many as possible weight matrix pairs satisfying Hurwitz-Radon (HR) orthogonality. 
We would, of course, like all the weight matrices to satisfy HR-orthogonality, but there is a limit to 
this number [1] which, except for the Alamouti code, is much lesser than 2Tnr, the number of weight 
matrices of a full-rate STBC when rir < nt. It can easily be checked that for the Alamouti code, I2 = C2. 
It is known that for a rate-1 code for nt > 2, one cannot have all the full-ranked weight matrices mutually 
satisfying HR-orthogonality. For such STBCs, the minimum self-interference is achieved if the STBCs 
are (7-group decodable, with g as large as possible. At present, the best known rate-1 low complexity 
multi-group decodable codes are the 4-group decodable codes for any number of transmit antennas [27], 
[28], [29]. These codes are not full-rate for n^. > 1. If one were to require a full-rate code, the codes 
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Ai = I„ 


A±+i 




A(g-l)fc 1 ^ 

1 


As 


"A ^ -i-2 — -As-A fc 1 




A(g-l)fc 1 r, = A2A(g-l)fc | -j^ 

a a 










Afc 

9 


A2fc = Afc Afc_|_-|^ 

9 9 




Ak = Afc A (g-l)fc 1 
9 g ^ 



lABUi i 

Weight matrices of a o-group decodable code 



in literature [27], [28], [29] are not suitable for extension to higher number of receive antennas, since 
their design is obtained by iterative methods. In the next section, we propose a new design methodology 
to obtain the weight matrices of a rate-1, 4-group decodable code by algebraic methods for 2" transmit 
antennas. These codes can be extended to higher number of receive antennas to obtain full-rate STBCs 
with lower ML-decoding complexity and lower self-interference than the existing designs. 

IV. Construction of Rate-1, 4-group decodable codes 

We make use of the following theorem, presented in [28], to construct rate-1, 4-group decodable codes 
for n = 2" transmit antennas. 

Theorem 1: [28] An n x n linear dispersion code transmitting k real symbols is ^r-group decodable if 
the weight matrices satisfy the following conditions: 

1) Af = I„, iG{l,2,--- ,|}. 

2) A2 = -I„, j e{Tf + l,m = l,2,--- ,5-1}- 

3) A,A, =A,A„ i,iG{l,2,--- , |}. 

4) A,A, =A,Ai, {1,2,--- j g + 1, m = 1, 2, • • • , g - 1}. 

5) AiA,- = -A,A,, i,jE{2|=+l,m = l,2,-- - ,5-1}, i^j. 

6) A:^+. = AiA™fc+i, m G {1, 2, • • • , 5 - 1}, i G {1, 2, • • • , f }. 

99 ^ 

Table U illustrates the weight matrices of a (7-group decodable code which satisfy the above conditions. 
The weight matrices in each column belong to the same group. 

In order to obtain a rate-1, 4-group decodable STBC for 2"- transmit antennas, it is sufficient if we have 
2a+i matrices satisfying the conditions in Theorem [T] To obtain theseQ we make use of the following 
lemmas. 

^These STBCs can be obtained elegantly using the theory of Clifford Algebra but to make the paper accessible to a wider 
group of readers, we have preferred to make use of simple concepts from matrix theory without reference to Clifford Algebra. 
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Lemma I: [33] If n = 2° and invertible complex matrices of size nxn, denoted by Fj, z = 1, 2, • • • , 2a, 
anticommute pairwise, then the set of products Fj^Fj^ • • • Fj^ with 1 < ii < • • • < < 2a along with I„ 
forms a basis for the 2^°^ dimensional space of all n x n matrices over C. 

Proof: The proof is provided for the sake of completeness. Assume that in the set of products 
FjjFjj • • • Fj^, 1 < ii < • • • < is < 2a, along with I„, at most k elements are linearly independent over 
C, for some /c < 2^". So, 

k+l 

J]aiF^F2- •••F2-» =0„, a^^O, A,^. G {0, 1}, j = 1, 2, • • • , 2a. (6) 

i=l 

Noting that F2 • • • F2a anticommutes with Fi but commutes with each of F2, • • • , F2a, premultiplying 
each term of Q by F2 • • • F2a results in a new equation with the coefficients a-i negated for those terms 
in ^ containing Fi . Adding this new equation to Q yields another equation containing fewer summands 
than leading to a contradiction. So, k = 2^°^, which proves the theorem. ■ 
Lemma 2: If all the mutually anticommuting nxn matrices Fj,i = 1,2,- •• ,2a, are unitary and 
anti-Hermitian, so that they square to —In, then the product Fj^Fj^ • • • Fj^ with 1 < ii < • • • < < 2a 
squares to (—1) 2 I^. 
Proof: We have 

(F,,F,, • • •F,J(F,,F,, • • -F.J = (-1)^-1(f2 F,, • ■■FJiF.^F,, • • • F, J 

= (-1)-1(-1)-2(f2 p2 . . . F,J(F,3F,, • • • F,J 

= (_l)[(^-l)+(s-2)+-l]^p2 p2 ...^2^ 
= (-1)^(-1)^I„ = (-1)^I„, 

which proves the lemma. ■ 
Lemma 3: Let Fi,i = 1,2, ••• ,2a be anticommuting, anti-Hermitian, unitary matrices. Let Jli = 
{Fi^,Fi^,- ■ ■ ,FiJ and ^2 = {Fj,,Fj^, ■ ■ ■ ,FjJ with 1 < ii < • • • < < 2a and 1 < ji < • • • < < 
2a. Let \Qi n 0.2] = P- Then the product matrix Fj^Fj^ • • - Fj^ commutes with Fj^Fj^ • • - Fj^ if exactly 
one of the following is satisfied, and anticommutes otherwise. 

1) r,s and p are all odd. 

2) The product rs is even and p is even (including 0). 
Proof: When F^^ G J^i n il.2, we note that 

(Fiji, ■ ■■FjFj, = {-l^'F,,(F^J^. ■ ■■FJ 
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and when F^-, ^ n ^2, we have {Fiji^ ■ • -F^jFj, = (-l)^Fj, (F^^Fi, • • -F^J. Now, 

(F,,F,, • • • F,J(F,,F,-, • • • F,J = (-l)f(-i)(-l)(-rt«(F,,F,, • • • F,v)(F,,F,, • • • F, J 

= (-ir-P(F,,F,vF,v)(F.,F,,...F,J. 

Case 1) Since r, s and p are all odd, (^—ly^-v - 1. 

Case 2) The product rs is even and p is even (including 0). Hence (— = 1. ■ 
From Theorem [U to get a rate-1, 4-group decodable STBC, we need 3 pairwise anticommuting, anti- 
Hermitian matrices which commute with a group of 2"^^ Hermitian, pairwise commuting matrices. Once 
these are identified, the other weight matrices can be easily obtained. From [2], one can obtain 2a 
pairwise anticommuting, anti-Hermitian matrices and the method to obtain these is presented here for 
completeness. Let 



Pi = 



1 




J 




1 




,P2 = 


,P3 = 




-1 




_ i _ 




-1 



and A®"* = A(8)A{X)A---(8)A. The 2a anti-Hermitian, pairwise anti-commuting matrices are 

m times 

Fi = ±jPf", 

F2fc = lf"'(g)Pi(g)Pf"\ A: = l,...,a, 

F2fc+i = lf-^(g)P2(g)Pr\ fc = l,---,a-l. 

Henceforth, Fj, i = 1, 2, • • • , 2a, refer to the matrices obtained using the above method. For a set 5 = 
{ai,a2,--- ,an}, define F(5) as 

P(cS)^{a^a^^---a^,A, G{0,1}}. 

We choose Fi, F2 and F3 to be the three pairwise anticommuting, anti-Hermitian matrices (to be placed 
in the top row along with I„ in TableU Consider the set S = {jF4F5, JF6F7, • • • , jF2a-2F2a-i5 F1F2F3}, 
the cardinality of which is a — 1. Using Lemma |2] and Lemma [3l one can note that S consists of pairwise 
commuting matrices which are Hermitian. Moreover, it is clear that each of the matrices in the set 
also commutes with Fi, F2 and F3. Hence, P(5), which has cardinality 2"~^ is also a set with pairwise 
commuting, Hermitian matrices which also commute with Fi, F2 and F3. The linear independence of P(5) 
over M is easy to see by applying Lemma [1] Hence, we have 3 pairwise anticommuting, anti-Hermitian 
matrices which commute with a group of 2"~^ Hermitian, pairwise commuting matrices. Having obtained 
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Is 


Fi 


F2 


F3 


jF4F5 


iFiF4F5 


JF2F4F5 


jF3F4F5 


F1F2F3 


— F2F3 


F1F3 


— F1F2 


JF1F2F3F4F5 


-JF2F3F4F5 


JF1F3F4F5 


-iFiF2F4F5 



TABLE II 

Weight matrices of a rate-1, 4-group decodable STBC for 8 transmit antennas 



these, the other weight matrices are obtained from Theorem [T] To illustrate with an example, we consider 
n = 8 and show below how the weight matrices are obtained for the rate-1, 4-group decodable code. 

A. An example - n = 8 

Let Fj, i = 1, 2, • • ■ , 6 denote the 6 pairwise anticommuting, anti-Hermitian matrices. Choose Fi, F2 
and F3 to be the three anticommuting matrices required for code construction. Let 

S = {iF4F5,FiF2F3}, P(cS) = {l8,jF4F5,FiF2F3,jFiF2F3F4F5}. 

The 16 weight matrices of the rate-1, 4-group decodable code for 8 antennas are as shown in Table 
lU Each column corresponds to the weight matrices belonging to the same group. Note that the product 
of any two matrices in the first group is some other matrix in the same group. 

B. Coding gain calculations 

Let A(S,S') = (iet(ASAS^), where AS = S — S',S / S' denotes the codeword difference matrix. 
Let Asi = Si — s'^,i = 1,2, ■■ ■ , 2nt, where Sj and s[ are the real symbols encoding codeword matrices 
S and S', respectively. Hence, 

(2nt 2nt \ / 2nt 2nt \ 

i=l m=l / \i=l m=l / 

Note that because of the nature of construction of the weight matrices, we have 

AiA^ = AE^+iA^^^, |l,2,... ,^}, pG{l,2,3}. 

Further, since the code is 4-group decodable, 

I I ^ 2 2 

A(S,SO = det Yl Y Yl ^^i^s„^AiA 



H 

m 

p=0 \i=£^+l i=P^+l m=i+l 
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All the weight matrices in the first group are Hermitian and pairwise commuting and the product of any 
two such matrices is some other matrix in the same group. It is well known that commuting matrices 
are simultaneously diagonalizable. Hence, 



A, = ED,E^, i G |2,3,--- 



where Dj is a diagonal matrix. Since Aj is Hermitian as well as unitary, the diagonal elements of Dj are 
±1. The following lemma proves that Aj is traceless. 

Lemma 4: Let Fj,i = 1,2, ••• ,2a be 2'^ x 2°- unitary, pairwise anticommuting matrices. Then, the 
product matrix F^^F2^ • • • F2^", Aj G {0, 1}, i = 1, 2, • • • , 2a, with the exception of I20, is traceless. 

Proof: It is well known that tr(AB) = ir(BA) for any two matrices A and B. Let A and B be two 
invertible, n x n anticommuting matrices. Then, ABA"^ = — B. So, 

tr(ABA"^) = -ir(B) <^ tr{B) = -tr{B). 

.■.tr(B) = 0. (7) 

Similarly, it can be shown that tr{A) = 0. By applying Lemma [3l it can be seen that any product 
matrix F^^Fg^ • • • Fg^", excluding 12", anticommutes with some other invertible product matrix from the 
set {F^^Fj' • • • F2^°, Ai G {0, = 1, 2, 3, • • • , 2a}. Hence, from ©, we can say that every product 
matrix F^'F^" ■ ■ ■ Fg^" except is traceless. ■ 
From the above lemma, Aj except identity is traceless. Hence, Dj has an equal number of 'I's and '-I's. 
In fact, because of the nature of construction of the matrices Fi,i = 1, 2, • • • , 2a, the product matrices 
FjFj_|_i, for even i, and the product matrix F1F2F3 are always diagonal (easily seen from the definition of 
Fj, i = 1, 2, • • • , 2a). Hence, all the weight matrices of the first group excluding Ai = I„j are diagonal, 
with the diagonal elements being ±1. Since these diagonal matrices also commute with F2 and F3, the 
diagonal entries are such that for every odd i, if the {i, if^ entry is 1(-1), then, the (i + 1, i + 1)*'^ entry 
is also 1(-1, respectively). To summarize, the properties of Aj, z = 2, • • • , ^ are listed below. 

A._A^^. A^-I 

Ai(m,n) = 0, m / n; Ai(j,j) = ±1, j = 1, 2, • • • , nt, 

tr(A,) = 0, (8) 
A,(j,j) = A,(j + l,i + l), i = l,3,5,--- (9) 
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AjAj = Afc, G |l,2, ••• 



(10) 



In view of these properties, 



{P + I)it 1 (P + l)nt 



A(S,S') = deti^i E I I ' 

\p=0 \i=E^+l i=S^+l m=i+l 

where Dj^ = AjA^ = A^ for some G {1, 2, • • • , and 

rit 3 

j=lp=0 \i=l 



(11) 



where dij = ±1 and dij = 1. In fact, dij = Aj(j, j), i = 1, 2, 3, • • • , From ([TT]) . A(S, S') is a product 
of the sum of squares and it is minimized when only one group, say p = 0, gives a non-zero contribution. 
Hence, 



nt 



min(A(S,S')) = min JJ ^dij^Si 



7 = 1 \i=l 

where min(y) denotes the minimum value of y over all possible values of x. From 

X 



minf A(S, S')) = niin 

S,S' ^ ^ ' As, 



(12) 



^j=i \i=i 

We need the minimum determinant to be as high a non-zero number as possible. In this regard, let 



W 



y^Kj], Wij = di(^2j~l), = 1,2,--- ,y 



(13) 



and 



— [yitp+i ' yntE+2'' ' ' ' ' y ntAp+i) \ ~ ^['^itp+i i ' ■ ■ ■ ' ^ >»t(p+i) ] ' p — o, i, 2, 3. 



Lemma 5: W as defined in (fT3] ) is an orthogonal matrix. 

Proof: From ([T3] ). it can be noted that the columns of W are obtained from the diagonal elements 
of Aj, i = 1, 2, • • • , ^. Each element of a column i of W corresponds to every odd numbered diagonal 
element of Aj. Denote the i*^ column of W by W/. Applying ( fTOl ) and ([8]l in that order, 

1 1 
(wi,Wj) = — ir(AiAj) = — tr(Afc) = 
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where 



Hence, W is orthogonal. 

Substituting in ([T2l ). we get 



0, if i ^ j, 

1, otherwise. 



min(A(S,S')) = min 



yo 



So, the minimum determinant is a power of the minimum product distance in nt/2 real dimensions. If 
£ , the product distance can be maximized by premultiplying y^ with a suitable orthogonal rotation 
matrix V given in [34]. This operation maximizes the minimum determinant and hence the coding gain. 
So, the 2nt real symbols of the rate-1, 4-group decodable code are encoded by grouping ^ real symbols 
each into 4 groups and each group of symbols taking value from a unitarily rotated vector belonging to 
, the rotation matrix being W"^V. For 4 transmit antennas. 



W 



and for 8 transmit antennas. 



1 

V2 



1 -1 




0.8507 


-0.5257 




, v = 






1 1 




0.5257 


0.8507 



w 



-0.3664 


-0.7677 


0.4231 


0.3121 


-0.2264 


-0.4745 


-0.6846 


-0.5050 


-0.4745 


0.2264 


-0.5050 


0.6846 


-0.7677 


0.3664 


0.3121 


-0.4231 



If the practically used square QAM constellation of size M is used, encoding is done as follows : the nt 
complex symbols in each codeword matrix take values from the M-QAM and are split into two groups, 
one group consisting of the real parts of the nt symbols and the other group consisting of the imaginary 
parts. Each group is further divided into two subgroups, each consisting of nt/2 real symbols. So, in 
all, there are 4 groups consisting of nt/2 real symbols. As used before, denoting the column vectors 
consisting of the symbols in a group by y^, p = 0, 1,2,3 (the entries of y^ take values independently 
from \/M-PAM), let Sp = W^Vy^, where W and V are as explained before. Then the codeword matrix 
is given by 



EE' 

p=0 i=l 
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Note that the above codeword matrix can also be expressed as 

3 T" 

s = EEy^.AW+.' (14) 

p=0 i=l 

where Apr.t ^ . = V^-^i ^jjAp^t , j. p = 0, 1,2,3, with ujji being the element of W"^V. Clearly, 

the weight matrices A'pr>t p = 0, 1, 2, 3, satisfy the condition 

A'^^^ (a'^,,^.)"" + A'^^^. (Ak+,)'' = 0„„ for < / < p < 3, and j = 0, 1, • • • , f . 

nt -2 

Consequently, the ML-decoding complexity of the code is of the order of M 4 . This is because there 
are four groups consisting of nj/2 real symbols each and the symbols in each group can be decoded 
independently from the symbols in the other groups. In decoding the symbols in the same group jointly, 
one needs to make a search over \fM ^ = M~t possibilities for the symbols, since the real and the 
imaginary parts of a signal point in a square M-QAM have only ^/M possible values each (the real and 
the imaginary parts of a signal point of a square M-QAM take values from a \/M-PAM constellation). 
However, one need not make an exhaustive search over all the possible values for the nt/2 symbols. 
For every possible value of the first ^ — 1 real symbols, the last symbol is evaluated by quantization 

I ^-1 Tit-2 

[4]. Hence, the worst case ML-decoding complexity is of the order of vAf ^ = M 4 only. Fig. 
[J gives a comparison of the symbol error rate for the proposed STBC, the 4-group decodable STBC 
proposed by Yuen et al. [27] and the 4-group-decodable STBC proposed by Rajan [28], all for the 8x1 
MIMO system. The plots reveal that all the STBCs have the same performance for QAM constellations. 
Independently, we have computed that all the three codes have the same minimum determinant for QAM 
constellations. 

V. Extension to higher number of receive antennas 

When = 1, a rate-1, 4-group decodable STBC is the best full-rate STBC possible in terms of 
ML-decoding complexity and as a result, ergodic capacity. However, when > 1, we need more weight 
matrices to meet the full-rate criterion. In literature, there does not exist a 4-group decodable STBC with 
rate greater than 1. So, it is unlikely, though not proven, that there exists a full-rate, multi-group ML- 
decodable STBC with full-diversity for > 1. So, for > 1, we relax the requirement of multi-group 
decodability and simply aim for some reduction in the ML-decoding complexity and self-interference. 
Let Ui = 2°-. We know that if Fj, i = 1, 2, • • • , 2a are pairwise anticommuting, invertible matrices, then, 
the set T = {^i'^2' ' ' ' ^2a^ with G {0, 1}, i = 1, 2, • • • , 2a} is Unearly independent over C. Hence, 
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the set A4 = {T, jF} is linearly independent over M. As a result, the elements of M. can be used as 
weight matrices of a full-rate STBC for Ur > 1. Keeping in view that the self-interference has to be 
minimized, it is important to choose the weight matrices judiciously. The idea is that given a full-rate 
STBC for Ur — 1 receive antennas, obtain the additional weight matrices of a full-rate STBC for 
receive antennas by using the weight matrices of a rate-1, 4-group decodable STBC such that after the 
addition of the new weight matrices, the set consisting of the weight matrices of the rate-n^ code is 
linearly independent over M. This is achieved as follows. 

1) Obtain a rate-1, 4-group decodable STBC by using the construction detailed in Section [TVl Due to 
the nature of the construction, the product of any two weight matrices is always some other weight 
matrix of the code, up to negation. Denote the set of weight matrices by Qi. 

2) From the set T, choose a matrix that does not belong to Qi and multiply it with the elements of 
Qi to obtain a new set of weight matrices, denoted by C/2- Clearly, the two sets will not have any 
matrix in common. To see this, let A G and B G Tf]{A4 \ Qi), where B is the matrix chosen 
to be multiplied with the elements of Qi. Let BA = C G t/i. Hence, B = CA^ = itCA and CA 
belongs to Qi, up to negation. This contradicts the fact that B G !Ff\{M\Gi)- So, C cannot belong 
to Gi. 

The weight matrices of Q2 form a new, rate-1, 4-group decodable STBC. This is because the ML- 
decoding complexity does not change by multiplying the weight matrices of a code with a unitary 
matrix. In this case, we have multiplied the elements of ^1 with an element of T, which is a unitary 
matrix. Now, Qi [j Q2 is the set of weight matrices of a rate-2 code with an ML-decoding complexity 
of M"'.M 4 = M i . This is achieved by decoding the last nt symbols with a complexity of 
M"* and then conditionally decoding the first nt symbols using the 4-group decodability property 
as explained in Section IIV-BI 

3) For increasing n^, repeat as in the second step, obtaining new rate-1, 4-group decodable codes 
and then appending their weight matrices to obtain a new, rate-n^ code with an ML-decoding 
complexity of Af"*(""~4)^0-^. The new set of weight matrices is ljr=i 

4) When all the elements of T have been exhausted (this occurs when = nt/2), Step 3 can be 
continued till = nt by choosing the matrices that are to be multiplied with the elements of Qi 
from jTf]{A4 \ [j^li^ Gi)- Note from Lemma [T] that this does not spoil the linear independence 
of the weight matrices over R. 
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Note : In the case of the Perfect codes for nt transmit antennas, a layer [25], [26] corresponds to 
nt complex symbols. In case of our generalized Silver codes, a layer corresponds to a rate-1, 4-group 
decodable code encoding nt complex symbols. Also, the Silver code for an rit x system refers to 
the STBC containing rimin = min{7it,nr) individual rate-1, 4-group decodable codes, a property due to 
which self-interference is greatly reduced compared with other known full-rate codes. 

A. An illustration for nt = A 

For nt = 4, let Fi,F2,F3 and F4 be the four anticommuting, anti-Hermitian matrices obtained by 
the method presented in [2]. Let T = {Y^^Y^^'Y^^'YI^ \i G {0,1}, i = 1,2,3,4}. The rate-1, 4-group 
decodable code has the following 8 weight matrices, with weight matrices in each column belonging to 
the same group: 



I4 


Fi 


F2 


F3 


F1F2F3 


— F2F3 


F1F3 


— F1F2 



Hence, Qi = {I4, Fi, F2, F3, F1F2F3, — F2F3,FiF3, — F1F2}, Now, we choose a matrix from F which 
does not belong to ^1. One such matrix is F4. Pre-multiplying all the elements of Qi with Fi and applying 
the anticommuting property, we obtain a new rate- 1 , 4-group decodable code, whose weight matrices are 
as follows: 



F4 


— F1F4 


— F2F4 


— F3F4 


— F1F2F3F4 


— F2F3F4 


F1F3F4 


— F1F2F4 



Hence, ^2 = ^iQi = {F4,-FiF4, -F2F4, -F3F4, -F1F2F3F4, -F2F3F4, F1F3F4, -F1F2F4} and Gi 
[J Q2 is the set of weight matrices of the rate-2 STBC, which is full rate with an ML-decoding complexity 
of the order of M"^-^. 

Now, since there are no more elements left in F (neglecting negation), we can choose elements from 
jF. To construct a rate-3 code for 3 transmit antennas, we multiply the elements of Qi by JI4 to obtain 
the set ^3 = jQi. The weight matrices of the rate-3 code constitute the set Qi [J ^2 U^3- Similarly, the 
weight matrices of a full-rate code for > 4 are the elements of the set t/i IJ t/2 U ^3 U ^4' where 
^4 = jF^Qi = jC/2- It is obvious that Qi, Q2, G3 and Q4 represent the weight matrices of four individual 
rate-1, 4-group decodable codes, respectively. 
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B. Structure of the R-matrix and ML-decoding complexity 

The popular sphere decoding [36] technique is used to perform the ML-decoding of linear dispersion 
STBCs utilizing lattice constellations. A QR-decomposition of Hgg, the equivalent channel matrix, is 
performed to obtain H^q = QR and the ML-decoding metric is given by 



M(s) 



vec(Y) 



SNR 



HegS 



SNR 



Rs 



where y' = Q"^f ec(Y). The R-matrix of the Silver code for the nt x system has the following structure, 
irrespective of the channel realization: 



R 



D 



X 
D 



O 



2nt 



o 



2nt 



X 
X 

D 



where X G ]]j2ntx2nt ^ random non-sparse matrix whose entries depend on the channel coefficients 
and D = I4 (8)T, with T G M'a'^^ being an upper triangular matrix. The reason for this structure is that 
the weight matrices of the Silver code for an nt x system are also the weight matrices of min{nt,nj.) 
separate rate-1, 4-group decodable codes (as illustrated in Sec. |Vl). As a result of the structure of D, the 
R-matrix has a large number of zeros in the upper block, and hence, compared to other existing codes, 
the generalized Silver codes have lower average ML-decoding complexity. The worst case ML-decoding 
complexity is of the order of (M"'(""""~^))(M~^) = which is because in decoding 

the symbols, a search is to be made over all possible values of the last nt{nmin — 1) complex symbols 
(which requires a complexity of the order of M"' ~^)), while the remaining nt symbols can be 
conditionally decoded with a complexity of only, once the last ntyrimin — 1) symbols are fixed 

(a detailed explanation on conditional ML-decoding has been presented in [11], [4]). In simple words, to 
decode the Silver code, one does not need a 2njnmm dimensional real sphere decoder. All one requires is 
a 2nt{nmin — 1) dimensional real sphere decoder in conjunction with four parallel {nt — 2) /2 dimensional 
real sphere decoders. The decrease in the ML-decoding complexity is evident from the decrease in the 
dimension of the real sphere decoder from 2ntnmin to 2nt{nmin — 1) + ^^^^^5^- 
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C. Information Losslessness for rij. > rit 

For rir > nt, the Silver code is information lossless because its normalized generator matrix (normal- 
ization is done to ensure an appropriate SNR at each receive antenna) is orthogonal. To see this, the 
generator matrix for > 7it is given as 

[7;ec(Ai) ?;ec(A2) ••• vec{A2n^J], 



/nt 

where Aj G M^,i = 1,2,- •• ,2n^, are the weight matrices obtained as mentioned in Sec. |Vl with 
M = where F = {F^F^^ • • • F^-, G {0, 1}, i = 1, 2, 3, • • • , 2a}. For i,j G {1, 2, • • • , 2n2}, 

we have 

{vec{Ai),vec{Aj)) = real (tr (Af Aj)) (15) 

= ±real(ir(AiAj)) (16) 

real(tr(I„J) if i = j 

real (tr(jl„j) if Aj = jAj 

ibreal (tr(Afc)) otherwise, where ± A^ G \ {Inj , j'lnt } 
= ntSij. (17) 

Equation (fT6l ) holds because Aj, i = 1, • • • , 2n1 are either Hermitian or anti-Hermitian, and (fTTl ) follows 
from Lemma |4l 

Lemma 6: Tracelessness of the self-interference matrix is equivalent to column orthogonaUty of the 
generator matrix. 

Proof Using the definition of the self-interference matrix S*"*, given in Definition |6l 

2fc-l 2k 2k-l 2k 

tr (S*"*) = E E ii^^^f + ^^^^)] = 2 E E 'i'^ (^^^1 (A^Af )] ) 

1=1 j>i i=l j>i 

2fe-l 2k ^ 

= 2^^SiSj{vec{Ai),vec{Aj)), (18) 

where ( fTSl ) follows from ( fTSl ). From ( fTSl ). it is clear that column orthogonality of the generator matrix is 
equivalent to tracelessness of the self-interference matrix. ■ 
Recall that the second criterion given to maximize I2 (given by Q) requires that Sij = AjAj^+AjA^, 
i 7^ j, be traceless. It is clear from Lemma |6] that for our STBCs, Sij is traceless for i / j. 
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D. The Silver code for two transmit antennas 

The Silver code [7], [8] for two antennas, which is well known for being a low complexity, full-rate, 
full-diversity STBC for rij. > 2, transmits 2 complex symbols per channel use. A codeword matrix of 
the Silver code is given as 

•S5 + jSG S7 + jS8 
-S7+jS8 S^-jSQ 

where 

1+j l + 2j 
-l + 2j l-j 

The codeword encodes 8 real symbols si, S2, • ' ' ^sg, each taking values independently from a regular 
\/M-PAM constellation. The first four weight matrices are that of the Alamouti code, given by 





1 







j 










1 







3 


Ai = 






, A2 = 




, A3 = 






, A4 = 









1 







-j 




-1 







J 






Note that the Alamouti code is 4-group decodable for 2 transmit antennas. The Silver code's next 4 
weight matrices are obtained by multiplying the first four weight matrices by j. To make the code 
achieve full-diversity with the highest possible coding gain, post-multiplication by U is performed. It can 
be checked that U = -^(Ai+A2+A3 + 2A4). Effectively, the last 4 weight matrices of the silver code are 
jAjU, i = 1, - ■ ■ ,4, which also form another rate-1, 4-group decodable code. The unitary matrix U is so 
cleverly chosen that in addition to providing full-diversity with a high coding gain, the generator matrix 
is orthogonal (which can be checked using (fT6l)). making the code information lossless for > 2. The 
Silver code compares very well with the well known Golden code in error performance, while offering 
lower ML-decoding complexity of the order of M^. 

E. Achievability of Full-diversity 

The following theorem, (Theorem I, [35]) guarantees that full-diversity is possible for the generalized 
Silver codes with the real symbols taking values from PAM constellations, denoted by Apam- 

Theorem 2: For any given n x n square linear design 5 = |s = Yli=i ^i^i I ^ Apam, ^ = 1,2, 
• • • , k}, encoding k real symbols with full-rank weight matrices Aj, there exist G C, i = 1, • • • , A:, 
such that the STBC 5' = |s = X]i=i SjO^iAj | Sj G Apam, ^ = 1, 2, • • • , /cj offers full diversity. 

Since all the weight matrices of the generalized Silver code are either Hermitian or anti-Hermitian 
and hence full-ranked, full-diversity is achievable with the generalized Silver codes. However, finding 



Sl + jS2 S3 + jS4 
-S3+jS4, Si-jS2 



+ 3 



U 



1 

7f 
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out explicitly the values of Oj is an open problem. For the full-rate codes for 1 receive antenna, in 
Section IIV-BI we have identified the encoding scheme which not only provides full-diversity, but also 
maximizes the coding gain for PAM constellations. For the generalized Silver codes for higher number of 
receive antennas, each layer, corresponding to a rate-1, 4-group decodable code, is encoded as explained 
in IIV-BI Note from ([T4l) that this type of encoding neither reduces the number of matrix pairs satisfying 
Hurwitz-Radon orthogonality nor spoils the column orthogonality of the Generator matrix. In addition, 
we use a certain scaling factor to be multiplied with a certain subset of weight matrices to enhance the 
coding gain. The choice of the scaling factor is based on computer search. With the use of the scaling 
factor, the generalized Silver codes perform very well when compared with the punctured Perfect codes. 
Although we cannot mathematically prove that our codes have full-diversity with the constellation that 
we have used for simulation, the simulation plots seem to suggest that our codes have full-diversity, since 
the error performance of our codes matches that of the comparable punctured Perfect codes, which have 
been known to have full-diversity. 

VI. Simulation results 
In all the simulation scenarios in this section, we consider the Rayleigh block fading MIMO channel. 

A. 4 Tx 

We consider three MIMO systems - 4 x 2, 4 x 3 and 4x4 systems. The codes are constructed as 
illustrated in Subsection IV-AI To enhance the performance of our code for the 4x2 system, we have 
multiplied the weight matrices of Q2 (as defined in Subsection IV-Ab with the scalar e^'^^^. This is done 
primarily to enhance the coding gain, which was observed to be the highest when the scalar e^'^^^ was 
chosen. It is to be noted that this action does not alter the ML-decoding complexity and the column 
orthogonality of the generator matrix (so, the resultant weight matrices still satisfy the tracelessness 
criterion). Consequently, the weight matrices of the Silver code for the 4x2 system can be viewed to be 
from Gi U e^''''^g2. For the 4 x 3 MIMO system, the weight matrices of the Silver code are from the set 
Qi U e^'^/^Q2 y}jQi, while the weight matrices of the Silver code for the 4x4 system are from the set 
Qi y e^'^/'^Q2 yjjGi [J je^^^^G2- Fig- El shows the plot of the maximum mutual information achievable 
with our codes and the punctured Perfect codes [25] for 4 x 2 and 4x3 systems. In both the cases, 
our codes allow higher mutual information than the punctured Perfect code, as was expected. Regarding 
error performance, we have chosen 4 QAM for our simulations and encoding is done as explained in 
Subsection lEll 
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1) 4 X 2 MIMO 

Fig. [3] shows the plots of the symbol error rate (SER) as a function of the SNR at each receive 
antenna for five codes - the DjABBA code [12], the punctured Perfect code for 4 transmit antennas, 
the Silver code for the 4x2 system, the EAST code [37] and Oggier's code from crossed product 
Algebra with a provable NVD property [19]. Since the number of degrees of freedom of the channel 
is only 2, we use the Perfect code with 2 of its 4 layers punctured. Our code and the EAST code 
have the best performance. It is to be noted that the curves for the Silver code for the 4x2 system 
and the EAST code coincide. Also, the Silver code for the 4x2 system is the same as the one 
presented in [4], but has been designed using a new, systematic method. The Silver code for the 
4x2 system and the EAST code have an ML-decoding complexity of the order of M^'^ for square 
QAM constellation, while the DjABBA and Oggier's code have an ML-decoding complexity of 
order and M^-^, respectively. 

2) 4 X 3 MIMO 

Fig. m shows the plots of the SER as a function of the SNR at each receive antenna for two codes - 
the punctured perfect code (puncturing one of its 4 layers) and the Silver code for the 4x3 system. 
The Silver code for the 4x3 system has a marginally better performance than the punctured perfect 
code in the low to medium SNR range. It has an ML-decoding complexity of the order of A/®^ 
while that of the punctured Perfect code is M^^ (this reduction from M^^ to M^^ is due to the 
fact that the real and the imaginary parts of the last symbol can be evaluated by quantization, once 
the remaining symbols have been fixed). 

3) 4 X 4 MIMO 

Fig. H] shows the plots of the SER as a function of the SNR at each receive antenna for the Silver 
code for the 4x4 system and the Perfect code. The Silver code for the 4x4 system nearly 
matches the Perfect code in performance at low and medium SNR. More importantly, it has lower 
ML-decoding complexity of the order of M^^^, while that of the Perfect code is M^^. 

B. 8 Tx 

To construct the Silver code for the 8x2 system, we first construct a rate-1, 4-group decodable STBC 
as described in Section |IV] and denote the set of obtained weight matrices by Qi. Next we multiply the 
weight matrices of Qi by F4 to obtain a new set of weight matrices which is denoted by ^2- The weight 
matrices of the Silver code for the 8x2 system are obtained from Qi [j Q2- The Silver code for the 8x3 
system can be obtained by multiplying the matrices of Qi with Fg and appending the resulting weight 
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matrices to the set Qi[jQ2- The rival code is the punctured perfect code for 8 transmit antennas [26]. 
The maximum mutual information plots of the two codes are shown in Fig. |6] As expected, our code 
has higher mutual information, although lower than the ergodic capacity of the corresponding MIMO 
channels. 

Fig. |7] shows the symbol error performance of the Silver code for 8 x 2 system and the punctured 
Perfect code [26]. The constellation employed is 4-QAM. Again, to enhance performance by way of 
increasing the coding gain, we have multiplied the weight matrices of Q2 with the scalar e^, as done 
for the codes for 4 transmit antennas. The simulation plot suggests that our code has full diversity. The 
most important aspect of our code is that it has an ML-decoding complexity of M^'^, while that of the 
comparable punctured Perfect code is M^^. 

VII. Discussion 

In this paper, we proposed a scheme to obtain full-rate STBCs for 2" transmit antennas and any 
number of receive antennas with the lowest ML-decoding complexity and the least self-interference 
among known codes. The STBCs thus obtained allow higher mutual information than existing STBCs 
for the case Ur < fit. Identifying explicit constellations which can be mathematically proven to guarantee 
full-diversity and a non-vanishing determinant without increasing the ML-decoding complexity is an 
open problem. Also, one can seek to obtain full-rate STBCs with reduced ML-decoding complexity for 
arbitrary number of transmit (not a power of 2) and receive antennas. These are some of the directions 
for future research. 
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Fig. 1. SER comparison of the proposed STBC withi a few known 4-group decodable STBCs for tfie 8x1 MIMO system 
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Fig. 2. Ergodic capacity Vs SNR for codes for 4 x 2 and 4x3 systems 



April 4, 201 1 



DRAFT 



28 




10 ^' ' ' ' ' ' 

8 10 12 14 16 18 20 

SNR indb 



Fig. 3. SER performance at 4 BPCU for codes for 4 x 2 systems 
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Fig. 4. SER performance at 6 BPCU for codes for 4 x 3 systems 
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Fig. 5. SER performance at 8 BPCU for codes for 4 x 4 systems 
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Fig. 6. Ergodic capacity Vs SNR for codes for 8x2 and 8x3 systems 
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